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Let U m be an m x m Haar unitary matrix and /7r mn i be its n x n truncation. 
In this paper the large deviation is proven for the empirical eigenvalue density 
of C/r m>n ] as m/n — > A and n — > oo. The rate function and the limit distribution 
are given explicitly. U[ mjn \ is the random matrix model of quq, where u is a Haar 
unitary in a finite von Neumann algebra, q is a certain projection and they are 
free. The limit distribution coincides with the Brown measure of the operator 
quq. 
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1 Introduction 

b , 

Although the asymptotics of the eigenvalue density of different random matrices has been 
widely studied since the pioneering work of Wigner |18j . the first large deviation theorem 
for the empirical eigenvalue density of self-adjoint Gaussian random matrices was proven 
by Ben Arous and Guionnet much later [T]. After the publication of their work, several 
similar theorems were obtained for different kind of random matrices. In particular, Haar 
distributed unitaries were discussed by Hiai and Petz [TT] and the monograph [TU] contains 
more information about similar results (see also [T?^l2*]). Free probability theory has inspired 
non-commutative large deviation results for random matrices recently, see [Zj, for example. 

'•February 1, 2008 
2 E-mail: petz@math.bme.hu 
E-mail: refFyj@math.bme.hu 



1 



The aim of this article is to prove the large deviation theorem for the empirical eigenvalue 
density of truncated Haar unitary random matrices, and to determine the limit measure. 
Let U be an m x m Haar distributed unitary matrix. By truncating m — n bottom rows and 
m — n last columns, we get an n x n matrix. The truncated matrix is a contraction, hence 
the eigenvalues are in the unit disc. Our aim is to study the asymptotics of the empirical 
eigenvalue density when n — > oo and m/n — > A. The truncated Haar unitaries appeared 
in the works [T5J Ej. Since our random matrix model is unitarily invariant, the limiting 
eigenvalue density is rotation invariant in the complex plane. It turns out that the limiting 
density is supported on the disc of radius 1 / y/X. In the paper the large deviation result is 
established and the exact form of the rate function is given. The large deviation implies 
the weak convergence of the empirical eigenvalue density of the truncated unitaries with 
probability one. 

The paper is organized as follows. Section 2 contains some preliminaries about potential 
theory and large deviations. The large deviation result is stated in Section 3. Section 4 
contains the proof of our main result. In Section 5 of the paper we make a connection to 
free probability theory. The truncated Haar unitaries form a random matrix model for the 
non-commutative random variable quq, where u is an appropriate unitary, q is a projection 
and they are assumed to be free. We observe that the limiting eigenvalue density coincides 
with the Brown measure of the operator quq. Our paper is based on the joint eigenvalue 
density of truncated unitaries. In the Appendix we sketch the derivation of this formula 
following the original paper [19]. 

2 Preliminaries 

In this section we review the setting of large deviation for the empirical eigenvalue density 
of random matrices and collect some useful concepts and results from potential theory. 

Assume that T n (uS) is a random n x n matrix with complex eigenvalues Ci(u;), • • • > CnO^)- 
(If we want, we can fix an ordering of the eigenvalues, for example, regarding their absolute 
values and phases, but that is not necessary.) The empirical eigenvalue density of T n (u) is 
the random atomic measure 



where 5(z) denotes the Dirac measure supported on {z} C C. Therefore P n is a random 
measure, or a measure-valued random variable. 

Let us recall the definition of the large deviation principle [E] . Let (P n ) be a sequence of 
measures on a topological space X. The large deviation principle holds with rate function 
I : M(X) -> R+ U {+00} in the scale n~ 2 if 



P n (u) : 



tf(CiM) + . . . + 



n 




for all open set Gel, and 



limsup — logP n (F) < — inf I(x) 

n — >oo Tl x£F 
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for all closed set F G X. 

Let U(m) be an m x m Haar distributed unitary matrix. By truncating m — n bottom 
rows and m — n last columns, we get a n x n matrix £/[ m .„]. The truncated matrix U\ m ,n\ 
is not a unitary but its operator norm is at most 1. Hence the eigenvalues 0, (2, ■ ■ ■ , (n he 
in the disc T> := {z G C : \z\ < 1}. The relevant topological space is Ai(T>), the space of 
probability measures on T>. Note that this space is a compact metrizable space with respect 
to the weak convergence of measures. Let P\ m ,n\ be the empirical eigenvalue density of U[ m>n \. 
Hence P\ m ,n] ma Y be regarded as a measure on M.(V). 

We are going to benefit from the fact that the joint probability density of the eigenvalues 
of [/[ TOi „] is 

n 

7^- n lo-c/n^-iorr--- 1 

1 ' 1 l<i<j<n 1=1 



according to [12], see also the Appendix. The normalizing constant 

m — n + j — 1\ 1 



f= \ J J m-n + j 

was obtained in [Tlj . 

Next we recall some definitions and theorems of potential theory [Ej. For a signed measure 

S(z/) := y y log |z — it? I dv(z) dv{w) 
is the negative logarithmic energy of v. Since 

S(z/) = inf / / maxflog \z — w\, a) dviz) dviw), 
a<0 J Jv? 

this functional is upper semi-continuous. We want to show its concavity. 

The following lemma is strongly related to the properties of the logarithmic kernel K(z, w) = 
log \ z — w\ (cf. Theorem 1.16 in |T2]). 

Lemma 2.1 Let u be a compactly supported signed measure on C such that u(C) = 0. Then 
Tt(y) < ; and S(z/) = if and only if v — 0. 

From this lemma we can deduce strictly concavity of the functional E. First we prove that 



for all fJ-i, 1^2 ^ -M(P), moreover the equality holds if and only if ^ = /i 2 . For this, apply 
Lemma 12. II for the signed measure v = Hi — [12- The strict midpoint concavity (0) implies 
strict concavity by well-known arguments. 

Let K C C be a compact subset of the complex plane, and Ai(K) be the collection of all 
probability measures with support in K. The logarithmic energy E(fi) of a /i G M.(K) is 
defined as ^ 

log 1 7 dfi(z) dfi(w) 



w\ 



and the energy V of K by 

V:=w£{E(n) :|i£ M{K)}. 

The quantity 

cap(i^) := e~ v 

is called the logarithmic capacity of K. The logarithmic potential of/i G M {K) is the function 

U»:= f log . 1 . dy.{w) 



K 



\z — w\ 



defined on K. 



Let K C C be a closed set, and Q : X — > (— oo, oo] be a lower semi-continuous function 
which is finite on a set of positive capacity. The integral 



1 r 

log i r <i/x(,z) dfi(w) + 2 / Q{z)dji(z) 

K 2 \Z — W\ J K 



is called weighted energy. 

The following result tells about the minimizer of the weighted potential (cf. Theorem 1.3.3 



m 



[i5j; 



Proposition 2.2 Let Q as above. Assume that o G M.(K) has compact support, E{a) < oo 
and 

U a (z) + Q(z) 

coincides with a constant F on the support of a and is at least as large as F on K. Then a 
is the unique measure in M.(K) such that 

i q(°)= ^f J Q (/i), 
i.e., a is the so-called equilibrium measure associated with Q. 

The following lemma is the specialization of Proposition 12.21 to a radially symmetric func- 
tion Q : V — >• (—00,00], i. e., Q(z) = Q(\z\). We assume that Q is differentiable on (0,1) 
with absolute continuous derivative bounded below, moreover rQ'(r) increasing on (0, 1) and 

lim rQ'(r) = 00. 

r— >1 

Let ro > be the smallest number for which Q'(r) > for all r > vq, and we set Ro be the 
smallest solution of RqQ'(R ) = 1. Clearly < r < Rq < 1. 



Lemma 2.3 If the above conditions hold, them the functional Iq attains its minimum at a 
unique measure [Aq supported on the annulus 

S Q = {z : r < |^| < -R }, 

and the density of hq is given by 

dfi Q (z) = —(rQ'(r))'drd(p, z = re llf> . 
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Proof. The proof is similar to the one of Theorem IV.6.1 in [15J. Using the formula 

1 f 27T , 1 , f — logr, if \z\ < r 

— / log - — - oIlo = < I,, -r i 

2tt J \z-re^\ [ -\og\z\, if \z\ > r, 

we get that 

i fRo i 

= Q(R,)-logRo-Q(z), 
for z G Sq, since r = or Q'(r ) = 0. We have 

U"(z) + Q(z)=Q{R )-]ogR , 

which is clearly a constant. 

Next we check that U^(z) + Q(z) > Q(Ro) — log-Ro for \z\ < r and for \z\ > Rq. So fig 
satisfies conditions of Theorem 12.21 and it must be the unique minimizer. □ 

3 The large deviation theorem 

Our large deviation theorem for truncated Haar unitaries is the following. 

Theorem 3.1 Let U\ myn i be the n x n truncation of an m x m Haar unitary random matrix 
and let 1 < A < oo. If m/n — > A as n —>■ oo, then the sequence of empirical eigenvalue 
densities P n = P\ m n ] satisfies the large deviation principle in the scale 1/n 2 with rate function 

I{n) := — / / log \z — w\ d/i(z) d/i(w) — (A — 1) / log(l — \z\ 2 ) d/i(z) + B, 
J Jv 2 Jv 

for /i G M.{V), where 

A 2 log A A 2 log(A-l) log(A-l) A-l 

d := . 

2 2 2 2 

Furthermore, there exists a unique /i G M.{V) given by the density 

(A-l)r 

dHo(z) = — -2 dr dip, z = re llp 

7r (1 — r z ) 

on {z : \z\ < 1/y/X} such that I(no) = 0. 



Set i 

F(z, w) := — log \z — w\ — (log(l — \z\ 2 ) + log(l — \w\ 2 )) 

and 

F a (z,w) := mm(F(z, w), a), 
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for a > 0. Since F a (z,w) is bounded and continuous 

/i £ h J J F a (z,w) djj(z) d/j(w). 

is continuous in the weak* topology, when the support of fi is restricted to a compact set. 
The functional I is written as 



/(/i) = J J F(z, w) dfi(z) d[i(w) + B 

= sup / / F a (z,w) d[i(z) d[i(w) + B . 
«>o J Jv 2 



hence / is lower semi-continuous. 

We can write / in the form 

I {H) = - (A - 1) / log(l - \z\ 2 ) d/i(z) + B. 

Jv 

Here the first part —£(//) is strictly convex (as it was established in the previous section) 
and the second part is affine in /x. Therefore / is a strictly convex functional. 

If X is compact and A is a base for the topology, then the large deviation principle is 
equivalent to the following conditions (Theorem 4.1.18 in 

—I(x)= inf | limsup — logP n (G) 1 = inf < liminf — logP„(G) 



x&G,GeA [ n ->oc 71 J xGG,GeA [ n-»oo n z 

for all x G X. We apply this result in the case X = A4(T>), and we choose 



fjf G A4(X>) : 



z kl z k2 dfi'(z)- / z fcl ^ fc2 d/i(z) 
© Jv 



< e for fci + ki < m 



to be G(fi; m, e). For /i G .M(X>) the sets m, e) form a neighborhood base of \i for the 
weak* topology of A4(T>), where m G N and e > 0. To obtain the theorem, we have to prove 
that 

—I{yi) > inf I limsup — logP n (G) 

G { n— >oo 71 



-/(//) < inf (liminf -\ log P n (G)\ 

G [ n-s-oo n z J 



where G runs over neighborhoods of //. 

The large deviation theorem implies the almost sure weak convergence. 

Theorem 3.2 Let U\ m ,n], Pn and /iq as in Theorem \S.l\ Then 

weakly with probability 1. 

The proof is standard, one benefits from the compactness of the level sets of the rate 
function and the Borel-Cantelli lemma is used, see 
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4 Proof of the large deviation 

In this section we prove Theorem 13.11 Our method is based on the explicit form of the joint 
eigenvalue density. 

First we compute the limit of the normalizing constant C[ TOiTl ] given in (jTj). 
B =: lim -^logCL,-] 



n-^oo 7T, 2 



n-1 , . 

m — n + j — 1 



.7=1 V 



J 



1 ^-^ n — 1 — i m — n — 1 + z 

= — lim > log . 

n^oo n — 1 *-~< n — 1 i 

Here the limit of a Riemannian sum can be recognized and this gives an integral: 

B = -J (l-x)log ( X ~l + X ) dx 

A 2 log A A 2 log(A- 1) log(A - 1) A- 1 
2 + 2 2 + 2 

The lower and upper estimates are stated in the form of lemmas. 
Lemma 4.1 For every fi G Ai(T>), 

inf | lim sup — logP n (G) > < — / / F(z,w) dfi(z) dfi(w) — B 

G { n-+oo n J J Jt> 2 

where G runs over a neighborhood base of /i. 

This is the easier estimate, one can follow the proof of the earlier large deviation theorems, 
see [110 OH. 

Lemma 4.2 For every /x 6 A4(V), 

inf | lim inf — log PJG) \ > - [ [ F(z, w) dfi(z) d^{w) - B, 

G [ n-oo n Z J J J V 2 

where G runs over a neighborhood base of /x. 
Proof. If 

J J F(z,w) d/i(z) dfi(w) 

is infinite, then we have a trivial case. Therefore we may assume that this double integral is 
finite. 

Since F(z, w) is bounded from below, we have 

F(z,w) d/i(z) dfj,(w) — lim / / F(z,w) dfj,k(z) d^k{w) 
v 2 fc - > °° J Jv 2 



with the conditional measure 

ti(Bnv t ) 

for all Borel set B, where 

V k :- > : \z\ < !--!>. 



k 

So it suffices to assume, that the support of ji is contained in T>k for some k EN. 

Next by possible regularization of the measure /i, we we may assume that fi has a contin- 
uous density / on the unit disc V, and 5 < f(z) for some 5 > 0. 

We want to partition the disc into annuli of equal measure. Let k = [y/n\, and choose 

= rW<r{ B) <...<ra<rW = l, 



({z = re^ : r G [rg, rf n) ]}) = i for 1 < i < A:. 



such that 

// ( : = /-c' r : /• e [/•; 

Note that 

k 2 <n< k(k + 2), 

and there exists a sequence l±, . . . ,1^ such that fc < /j < k + 2, for 1 < i < k, and X)i=i h = n - 
Now we partition radially. For fixed % let 

= ^ <^<-..< < ri B) = 27T, 



such that 

/.({^re^retr'irf'Ue^Uf]})^ for 1 < j < h 



In this way we divided T> into n pieces, S!"',...,^. Here 

*(lZf 5 )<* = / ^<^<i^, (3, 

i 

where £ n = 2/(y / n + 2) — > and £^ = l/{-\/n — 1) — > as n — > oo. We can suppose, that 

lim ( max diam ( ) ) = 0. (4) 

n^oo Vl<i<n V / J 

In each part we take a smaller one -D^, similarly to by dividing the radial and 
phase intervals above into three equal parts, and selecting the middle ones, so that 

5(1 - e n ) 



9n 

We set 



J D M 9n5 



(5) 



A n :={(Ci,...,Cn):C<e D™,l<i<n} 
For any neighborhood G of /i 

A n E{CEV n :^E G} 
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for every n large enough. Then 

Pn(G) > V n (A n ) 



1 J . . . jf exp ((n - 1) |> - 1) log (l - |C| 2 ) j 



l<i<j'<n 



X 



n L jrvj^-^ 



Here for the first part we establish 



ton (n-l)(A-l) A m . n (1 _ 2) 



A - 1 

lim \ min log (l — \C\ 2 ) 



= (A-1) / log (1 - |C| 2 ) /(C) ^C, 
Jo 

because of (JH) and verify 

liminf— log [ min \( — rj\ 

l<i<j<ra 

> JJ^f(Of(v)^og\C- v \dCd v . (6) 

for the second part. □ 
The last step is to minimize I. Now we apply Lemma \2. 31 for 

Q(z) :=-^log(l-M 2 ) 

on T>. This function satisfies the conditions of the lemma. Hence the support of the limit 
measure /xo is the disc 

S\ = \z : \z\ < 



and the density is given by 

1 1 (A — l)r 

djiQ = —(rQ(r))drdtp = ^drdip, z = re 1!p . 

7r 7T (1 — r 2 ) 



For this /Jo again by [T 



A-l, M ^ 1, , (A - l) 2 f 27T /"7a r log(l - r 2 ) 

log(A-l) + — log A - V _ / / 2 v drd<p 



2 2A & 2tt Jo y (1-r 2 ) 

A 1 log(A - 1) + ^rlogA - ^Alog (^-] + lV £ = 0. 



2 ' 2A ° 2 V V A 

The uniqueness of /io satisfying I(/io) = follows from the strict convexity of /. 

5 Some connection to free probability 

Let Q m be an m x m projection matrix of rank n, and let U m be an m x m Haar unitary. 
Then the matrix Q m UmQm has the same non-zero eigenvalues as U[ m , n ], but it has m — n zero 
eigenvalues. The large deviation result for U[ m>n ] is easily modified to have the following. 

Theorem 5.1 Let 1 < A < oo and Q m ,U m as above. If m/n — > A as n —>■ oo, then the 
sequence of empirical eigenvalue densities Q m U m Q m satisfies the large deviation principle in 
the scale 1/n 2 with rate function 



Furthermore, the measure 



Kji), if ii = (i - a- 1 )^ + A-y 

+oo, otherwise 



IMi = (1 - A-^Jo + A~Vq 
is the unique minimizer of I, and I{jxq) = 0. 



Now let M. be a von Neumann algebra and r be a faithful normal trace on A4. The pair 
(Ai,r) is often called a non- commutative probability space. A unitary u G A4 is called a 
i/aar unitary if r(M fc ) = for every non-zero integer k. Let q G .M be a projection such that 
r(g) = A. If u and g are free (see [10 or 17 for more details about free probability), then 
the above (U m , Q m ) is random matrix model of the pair (u, q). This means that 

-E (TrV(U m , U* m , Q m ) ^ r (V(u, u*, q)) 
m 

for any polynomial V of three non-commuting indeterminants. This statement is a particular 
case of Voiculescu's fundamental result about asymptotic freeness ( [TO] , or Theorem 4.3.5 on 
p. 154 in pi 



For an element a of the von Neumann algebra A4, the Fuglede-Kadison determinant can 
be defined by: 

A(a) := lim o expr (log(a*a + e/) 1/2 ) . 
It was shown by L.G. Brown in 1983 that the function 

A i-> 7 ^-logA(a - XI) 
2n 
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is subharmonic and its Laplacian (taken in the distribution sense) is a probability measure 
fi a concentrated on the spectrum of a jl]. This measure is called the Brown measure and it 
is a sort of extension of the spectral multiplicity measure of normal operators: 

r (g( a )) = / g( z ) dfx a (z) (7) 

Jc 

for any function g on C that is analytic in a domain containing the spectrum of a. The 
Brown measure is computed for quite a few examples in the paper 

Let u be a Haar unitary, and q = q* = q 2 be free from u. Then uq is a so-called R-diagonal 
operator and its Brown measure is rotation invariant in the complex plane. According to [Hj 
the Brown measure has an atom of mass 1 - A" 1 at zero and the absolute continuous part 
has a density 

^ r 9 dr dip (z = reP) 
7rA(l-r 2 ) 2 Y V ; 

on {z : \z\ < 1/\^X}. We just observe that this measure coincides with the limiting measure 
in our large deviation theorem. In the moment we cannot deduce the Brown measure from 
the large deviation result but it is definitely worthwhile to study the relation. 



Appendix 



Let U m be an m x m Haar unitary matrix and write it in the block-matrix form 



A B 
C D 



where A is an n x n, B is n x (m — n) , C is (m — n)xn and D is an (m — n) x (m — n) matrix. 
The space of n x n (complex) matrices is easily identified with K 2n and the push forward 
of the usual Lebesgue measure is denoted by X n . It was obtained in [Hj that for m > 2n, the 
distribution measure of the n x n matrix A is absolute continuous with respect to X n and 
the density is 

m— In -i 



C(n } m) det(l - A*A) m - 2n l m < l d\ n (A) . 



To determine the joint distribution of the eigenvalues Ci, C2, • • • , Cn of A, we need only the 
matrices A and C, and by a unitary transformation we transform A to an upper triangular 
form 

Ai, 3 
A 2j3 



/ Ci 




Al, 2 

C2 









c 2 





c 3 



A lin \ 
A2,n 



Cn 



(9) 



/ 



where Ci, C2, ■ ■ ■ , C n are the column vectors of the matrix C. It is well-known that the 
Jacobian of this transformation is a multiple of 



n 10—01 
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Note that the columns of the matrix (jOJ are normalized and pairwise orthogonal. Following 
the idea of [TH], we integrate out the variables A^j, A 2) i, . . . , Aj.^, Cj, i < n. 

One can construct (n — m) x (n — m) matrices X^ such that 



We have X« = I and 



Since 



x«=/ + ]Tx«^|x( fc ). 

Isfcl 



(10) 



fc<i 



c;c i + y %2A ki A ki = c;x®c i 



the vectors Cj satisfy the equations 

c*x«a = i - loi 2 . (ii) 

Geometrically, the point (Cu, . . . , C m _ n)i ) lies in the ellipsoid given by To compute the 
volume of this ellipsoid it is enough to know the determinant of X^ and this is obtained 
from the above recursion: 

ls»-i| jKl 1^1 

After this preparation we move to integration. First we integrate with respect to the last 
column Ai >n , A 2 ,„, . . . , A n _i in C„. For fixed A 1>n . . . A n _ ljTl the distribution of C 1>n , C m _ n „i . 
is uniform on the set 

|Ci jn | 2 + . . . + |C m _ n _i jn | 2 < 1 — |Cn| 2 — |Ai jra | 2 . . . |A„_i jra | 2 , 

i.e. inside the ellipsoid defined by ([lljh The volume of this m — n — 1 dimensional complex 
ellipsoid is 

/i I/- \2\in— n— 1 

l 1 ISnj ) _ /-. _ I > |2\m-n-l TT I/- |2 

Hpt. Yin) ~ V \^\ ) 11 M > 



Integration out of Aj )Tl gives a factor |Ci|~ 2 from (|TU|) and all together we obtain a factor 
(1 — |Cn| 2 ) m_n_1 from the last column. The same procedure may be applied to the other 
columns. 
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